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Impurity Effect on Superconducting Properties in Molecular Substituted 
Organic Superconductor k-(ET) 2 Cu(NCS) 2 
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We report an impurity effect in the organic superconductor k-(ET)2Cu(NCS)2 by substi- 
tution of the ET molecule with an analogue, bis(methyleneditio)tetrathiafulvalene (MT). The 
superconducting transition temperature decreases with increasing substitution. The in-plane 
magnetic penetration depth is enhanced with substitution, which is quantitatively attributed 
to the decrease in the in-plane mean free path. The enhancement of the penetration depth can 
also explain the reduction of the effective pinning in terms of the condensation energy. 
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Organic superconductors k-(ET)2Y (ET denotes 
bis(ethylenedithio)tetrathiafulvalene, X = Cu(NCS)2 
and Cu[N(CN)2]Br) consist of alternating layers of con- 
ducting ET donors and insulating X anions. 1 Reflecting 
the layered structure, the system is well characterized 
by highly anisotropic (quasi-two-dimensional) electronic 
properties. A prominent feature of the organic super- 
conductors is its clean system with less disorder. In k- 
(ET) 2 Cu(NCS)2, one of the most vigorously investigated 
organics, the in-plane mean free path Z|| grows to the 
length of approximately 150 - 200 nm. 2,3 Taking account 
of its long in-plane penetration depth (A|| ~ 0.5 /mi) and 
short in-plane coherence length (£|| ~ 5 nm), the lo- 
cal clean limit approximation 4 (£|| <C L <C An) can be 
properly applied to the superconducting state. 

It has been known that strong electron correlations 
play an important role on the anomalous superconduct- 
ing properties, such as existence of an antiferromagnetic 
phase next to the superconductivity. 5 This is presum- 
ably reminiscent of unconventional superconductivity in 
high-Tc cuprates, in which a d-wave symmetry of the or- 
der parameter has been widely accepted. However in or- 
ganic superconductors, although many efforts have been 
achieved and its unconventional (d-wave) symmetry may 
be predominantly accepted, no clear consensus has been 
reached yet. 6 

To identify the symmetry of the order parameter, a 
very helpful method would be to investigate impurity ef- 
fects. There have been many studies to introduce impu- 
rity (or disorder) into organics: anion 7-10 or donor 11-15 
substitution, fast electron 16 or X-ray 17 irradiation, and 
fast-cooling. 18 As expected, they commonly suppress 
T c with increasing disorder regardless of its origin. A 
fully deuterated ET donor (dET) substituted salt k- 
[(ET)i_ x (dET) x ] 2 Cu(NCS) 2 is a good example of very 
weakly disordered systems, in which T c increases with x 
because of chemical pressure. 14, 15 Nevertheless, the scat- 
tering time barely takes minimum at around x — 0.5, in- 
dicating that the substitution gives a small scattering. 15 
For more detailed investigation, a systematic introduc- 



tion and control of stronger disorder would be indispens- 
able. 

In this paper, we present the substitution of a relative 
donor molecule bis(methylenedithio)tetrathiafulvalene 
(abbreviated as MT or BMDT-TTF, see lower inset of 
Fig. 1) for the ET molecule. In the substituted samples, 
T c monotonically decreases with increasing substitution. 
This substitution enhances in-plane magnetic penetra- 
tion depth, which is interpreted as an impurity effect on 
the basis of the clean local limit approximation. 

Single crystals of the MT substituted k- 
(ET) 2 Cu(NCS)2 were grown by a standard electro- 
chemical technique. We label the samples using the MT 
concentration cmt as, for example, "MT1%", which was 
prepared with 1 mg of MT and 99mg of ET dissolved into 
100 mL 1,1,2-trichloroethane (5% ethanol). Samples of 
MT0%, MT0.01%, MT0.1%, MTO.5%, MT1%, MT5%, 
MT10%, and MT15% were prepared. The crystalline 
shape and quality of the substituted samples were 
almost the same ones with the pristine MT0% sample. 
The magnetization measurements were performed using 
a SQUID magnetometer (Quantum Design, MPMS-XL). 
Each sample was cooled with rates of approximately 
— 100 K/min from room temperature (quenched) and 
—0.2 K/min from 100 K (slow-cooled). Temperature and 
magnetic field dependences of the magnetization were 
measured in zero field cooled condition. All the data 
for MT0% shown here were taken from the previous 
study. 19 

First of all, to quantify the amount of the actual 
substitution, we carried out FT-IR spectroscopy in k- 
(dET) 2 Cu(NCS) 2 at c MT = 5,10, and 15% (data 
not shown in the figure). A peak intensity of the 
molecular vibration mode in terminal ethylene groups 
has been used to determine the substitution ratio 
in K-[(ET)i_s(dETy 2 Cu[N(CN) 2 ]Br, 20 since the peak 
height of the vibration mode linearly alters with x. In 
the present case as well, the peak height in the pristine 
K-(dET) 2 Cu(NCS) 2 reduces to approximately 85±4% in 
MT15%, although no significant variation is observed at 
more dilute Cmt within experimental error. Taking into 
account the linear decrease in T„ with cmt as described 
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Fig. 1. Temperature dependence of the susceptibility of the MT 
substituted k-(ET) 2 Cu(NCS) 2 . The data of MT0% were taken 
from the previous work. 19 Filled and open symbols are obtained 
in the slow-cooled and quenched conditions, respectively. The 
upper inset shows T c as a function of the concentration of MT. 
The filled (open) diamonds denote re-(ET)2Cu(NCS)2 in slow- 
cooled (quenched) condition, and inversed triangles the slow- 
coolcd data in re-(dET)2Cu(NCS)2- Broken lines are guides for 
the eye. The lower inset displays ET and MT molecules. 
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Fig. 2. Irreversibility field as a function of T/T c in the MT substi- 
tuted k-(ET) 2 Cu(NCS) 2 in the slow-cooled condition for MT0% 
(circle), MT1% (triangle), and MT5% (square). The upper inset 
shows the magnetization curves at 2.0 K, and the lower one the 
magnetization hysteresis widths at T/T c = 0.25. Broken curves 
are guides for the eye. 



below, it is concluded that the actual value of the MT 
concentration is in the same order with cmt in the whole 
concentration prepared in the present work. 

Figure 1 shows the temperature dependence of the sus- 
ceptibility in k-(ET)2Cu(NCS)2 at several cmt- The data 
were taken in a magnetic field of 3 Oe perpendicular to 
the conduction plane. The demagnetization factor was 
corrected by means of an ellipsoidal approximation. The 
slow-cooled condition (filled symbols) gives almost a full 
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Fig. 3. Semi-logarithmic plots of the reversible magnetization 
curves as a function of magnetic field in the MT substituted 
k-(ET) 2 Cu(NCS)2 (a) in MT1% and (b) MT5%, respectively. 
Data were taken in the slow-cooled condition. Data points be- 
low iT; rr were omitted for clarity. Solid lines indicate the linear 
regions used to evaluate Am . 
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Fig. 4. In-plane magnetic penetration depth of the MT substi- 
tuted k-(ET) 2 Cu(NCS) 2 at c M t = 0% (circles), 1% (triangles), 
and 5% (squares). The filled (open) symbols indicate the data in 
slow-cooled (quenched) condition. Solid curves are obtained on 
the basis of a d-wave model (see text), and a broken curve is an 
s-wave (local clean) model. 



Meissner volume fraction, whereas in the quenched case 
(open symbols) a slight suppression of superconductiv- 
ity (still above ~ 90%vol) is found in cmt > 1%- We 
do not find any trace of magnetic impurity for all the 
samples, implying that the present substitution induces 
non-magnetic disorder. T c is defined as the intercept 
of the extrapolation of the normal and superconduct- 
ing states. The superconducting transition in MT5% be- 
comes broader than that in MT0% and MT1%. Neverthe- 
less, as shown in the inset of Fig. 1, T c linearly decreases 
with increasing cmt in both cooling rates. The slow- 
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Fig. 5. In-plane magnetic penetration depth as a function of the 
mean free path in k-(ET) 2 X (X=Cu(NCS) 2 , Cu[N(CN) 2 ]Br). 
The data for the latter one are taken from the literature. 19 Small 
open circles are calculated from eq. (1) with fitting parameters of 
A L (0) = 420 and 523 nm for X =Cu(NCS) 2 and Cu[N(CN) 2 ]Br, 
respectively. Broken curves are guides for the eye. 



cooled data in the MT substituted K-(dET) 2 Cu(NCS) 2 
also have the same slope as in k-(ET) 2 Cu(NCS) 2 . 

Next we move on to the magnetization curve measure- 
ments. The upper inset of Fig. 2 shows the M(H) curves 
obtained at 2.0 K in the slow-cooled condition, where 
M is the volume magnetization. The irreversibility field 
-f/irr, denoted by arrows, decreases with increasing cmt- 
The hysteresis width of the magnetization (AM) below 
iJj rr at a reduced temperature T/T c = 0.25 is displayed 
in the lower inset of Fig. 2. The suppression of AM obvi- 
ously occurs in high Cmt- This behavior reflects its weak- 
ened effective vortex pinning. The reduced temperature 
(T/T c ) dependence of H m is shown in the main panel of 
Fig. 2. As described above, H m is suppressed with in- 
creasing cmt- Thus, the vortex liquid phase (reversible 
region in H > i?i rr ) expands. This also indicates that 
the effective pinning reduces with Cmt- The reduction of 
the effective pinning will be attributed to decreasing the 
condensation energy, as discussed further below. 

We turn to deriving the in-plane penetration depth 
from the magnetization curves. Figure 3 shows the semi- 
logarithmic plot of the reversible magnetization in the 
slow-cooled condition. An estimation of A|| is given 
from linear slopes as described below; for type-II su- 
perconductor in H c i < < H <C H C 2, the Lon- 
don model ' represents the volume magnetization as 
M = -(a0 o /327T 2 A|) ln(H c2 (3/H), where O is the mag- 
netic flux quantum, a(= 0.70) a correction factor for 
vortex-core contribution, 22 and [3 a constant of order 
unity. According to this relationship, steeper slope in the 
M vs. In H plot gives shorter Am. The solid lines shown 
in Fig. 3 exhibit linear regions fitted to this relation be- 
tween ffi rr and H c2 . As can be seen, the slope of the 
lines decreases with increasing temperature, so that A|| 
ordinarily increases with T. 

The temperature dependence of Am evaluated above is 



shown in Fig. 4. Although there is no cooling rate depen- 
dence of Am(T) within experimental accuracy, a marked 
enhancement of A|| with increasing cmt is recognized. 

In order to estimate A||(0), we apply an extrapolation 
to K by using a d-wave model: 23 ' 24 A||(T) = A , (0) [1 - 
at — bt 3 }~°- 5 7 where A||(0), a, and & are fitting parameters. 
The solid curves in Fig. 4 are obtained by fitting to the 
slow-cooled data (filled symbols) below T/T c < 0.7 with 
A||(0) = 430, 445, and 465 nm, a = 0.12, 0.69, and 0.76, 
and b = 1.32, 0.41, and 0.56 for MT0%, MT1%, and 
MT5%, respectively. The parameters a and b are in the 
same order of magnitude as the theoretically expected 
ones (a = 0.6478 and b = 0.2 76), 24 except for b in MT0%. 
On the other hand, as overlaid in Fig. 4, an s-wave (local 
clean) model with A||(0) = 435 nm (broken curve) also 
seems to be agreeable to the experimental data in MT0%. 
In the following, we adopt A||(0) obtained from the d- 
wave model fits, although we cannot assert that A||(T) 
truly reflects the unconventional symmetry within the 
present study. 

Before discussing the variation of Ay(0) by the MT 
substitution, we mention an intrinsic impurity effect on 
the penetration depth of the pristine samples. As in our 
previous report, 19 Ay(0) in k-(ET) 2 A quantitatively re- 
flects its cleanness; on the basis of the local clean BCS 
approximation, A|| (0) is described as 

A||(0) = A L (0)(l+£o/7||)°- 5 , (1) 

where Al(0)[= cm* /Airen s {Q)] and £o(= ahvp/kBT c ) are 
the London penetration depth and coherence length for 
an ideally pure sample, respectively, and c is the light 
velocity, to* the effective mass, n s (0) the carrier density, 
vf the Fermi velocity, and a — 0.18. This indicates that 
A|| (0) becomes long with decreasing lu. We have observed 
A||(0) = 430 nm in k-(ET) 2 Cu(NCS) 2 , which is almost 
the same with the calculated value of Al(0) = 410 nm. 
Taking account of eq. (1), this is due to its long mean free 
path, i|| ~ 150 nm. 2 On the other hand, Ay(0) = 570 nm 
obtained in K-(ET) 2 Cu[N(CN) 2 ]Br, which is relatively 
larger than Al(0) = 400 nm, is attributed to short In (~ 
40 nm), implying an intrinsically dirtier system. More- 
over, its fast-cooling process induces much more disorder, 
resulting in a further enhancement of A|| (0) with decreas- 
ing Z|| as displayed in Fig. 5 (filled diamonds). 19 This can 
also be interpreted on the basis of eq. (1) (small circles in 
Fig. 5). Thus, the local clean model adequately describes 
the absolute value of A||(0) in these organics. 

The MT molecules, most probably working as impu- 
rity, will give rise to scattering quasiparticles and de- 
creasing Consequently A||(0) is varied via Z||, whereas 
Al(0) would not be influenced by the MT substitution. 
Very recently, the scattering time T|| is systematically 
studied in the present system by means of the de Haas- 
van Alphen oscillation effect. 25 Since lu — Vft», one can 
obtain lu , only if vf is independent of cmt- In this exper- 
iment, T|| = 3 ps in MT0%, giving l\\ — 135 nm by apply- 
ing vf = 4.5 x 10 6 cm/s, varies into t\\ = 1.6 ps (so /|| = 
72 nm) in MT1% and 1 ps (45 nm) in MT5%, whereas 
m* does not alter at all. 25 As a function of X\\(0) is 
plotted in Fig. 5 (filled symbols). We note that the local 
clean limit condition holds for the present case as dis- 
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Table I. Superconducting parameters in the MT substituted k- 
(ET) 2 Cu(NCS) 2 in the slow-cooled condition. 
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played in Table I, where the in-plane coherence length at 
K, £||(0), is estimated from the form i^^ 1 = + 17 . 

Now, by using /|| obtained from T||, one can evalu- 
ate A||(0) according to eq. (1) with a fitting parameter 
Al(0) = 420 nm, which is shown in Fig. 5 (small open 
circles). As described in this figure, the calculated A||(0) 
is fairly consistent with the experimental data. Thus the 
absolute value of A||(0) is quantitatively explained from 
the view of the impurity effect in the same manner in 
fast-cooling of K -(ET) 2 Cu[N(CN) 2 ]Br. 

We finally discuss the effective vortex pinning in terms 
of the superconducting condensation energy. The distri- 
bution of impurity sites originating from the MT substi- 
tution is very likely to be uniform in the whole bulk. The 
size of them is presumably in a molecular scale (a few A) 
which is much smaller than the vortex core diameter, 



2i 



13 nm. Thus, the introduction of MT molecules 



does not contribute to vortex pinning, but curiously re- 
duces the effective pinning as seen in the suppression of 
AM with cmt- Such behavior is very similar to the case 
of fast-cooling in K-(ET) 2 Cu[N(CN) 2 ]Br, in which AM 
critically reduces with cooling rate. Although a disorder- 
domain induced pinning model has been proposed, 26 this 
is not applicable to the present case. As a reasonable ex- 
planation on the effective pinning reduced by both meth- 
ods, the MT substitution and fast-cooling, which may 
originate from a similar mechanism, we focus on the su- 
perconducting condensation energy. Namely, the pinning 
potential for a normal-superconducting boundary (U p ) is 
described as the condensation energy of the vortex core, 
which is given by U p = (l/2)/i ff c 2 ■ tt£ 2 = 0g/(16^ O 7rA|), 
where (1qH c is the thermodynamic critical field expressed 
as //o-ffc = </ ) o/(2v / 27rA||tf||). We note that U p depends on 
only A||, and in addition, smaller U p will cause weaker 
effective pinning. As shown in Table I, U p decreases 
with cmt, which is consistent with the present reduc- 
tion of the effective pinning. The similar calculation in 
K-(ET) 2 Cu[N(CN) 2 ]Br leads to U p = 2.1 and 1.4 x 10" 13 
J/m 3 in the slow-cooled and quenched conditions, respec- 
tively. This is the same trend with the MT substitution. 
For further quantitative discussion on the vortex pinning, 
however, the origin of vortex pinning should be clarified 
in the future. 

In conclusion, we have observed an impurity effect by 
molecular substitution. The substituted MT molecules 
clearly work as impurity giving rise to scattering quasi- 
particles, and so T c monotonically decreases. The en- 
hancement of A || (0) , which is adequately interpreted as 



the impurity effect on the basis of the local clean limit 
approximation, results in the decrease in the condensa- 
tion energy. It can also account for the reduction of the 
effective pinning. 
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